Perturbative analysis of coherent quantum ratchets in cold atom systems 
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We present a perturbative study of the response of cold atoms in an optical lattice to a weak 
time- and space-asymmetric periodic driving signal. In the noninteracting limit, and for a finite 
set of resonant frequencies, we show how a coherent, long lasting ratchet current results from the 
interference between first and second order processes. In those suitable three-level model 

can account for the entire dynamics, yielding surprisingly good agreement with numerically exact 
results for weak and moderately strong driving. 
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I. INTRODUCTION 

In recent years considerable attention has been di- 
rected at producing directed transport in driven systems, 
in which the driving does not have a net bias. This type 
of "ratchet effect" PI is important from the viewpoint of 
technology in controlling the passage of particles, ranging 
from electrons to nanospheres, and also to more unusual 
applications such as understanding the operation of bio- 
logical molecular motors. 

The most well-studied form of ratchet physics arises 
from the interplay between dissipation and the driving 
potential. The paradigmatic example is given by a Brow- 
nian particle in a periodic potential The system is 
driven from equilibrium by periodically varying the po- 
tential, and a ratchet current is produced when the rel- 
evant space and time symmetries of the driven system, 
which would otherwise forbid the formation of a directed 
current [3[, are broken. 

Perhaps surprisingly, however, dissipation is not a nec- 
essary requirement for the production of a ratchet, and 
even in strictly Hamiltonian systems a ratchet current 
can be produced jUljJ. Cold atoms loaded into optical 
lattices have emerged as excellent candidates to study 
such coherent ratchet effects, as the level of dissipation 
can be controlled rather precisely [Tl|, and indeed such 
systems can be arranged to be essentially dissipation- 
free. Initial investigations of ratchet-like behavior con- 
centrated on the "quantum kicked rotor" , which can be 
realized in experiment by subjecting a gas of ultracold 
atoms to a pulsed optical lattice potential. This has al- 
lowed the detailed experimental investigation of quantum 
chaos effects such as dynamical localization Q and quan- 
tum resonances 0] which can be harnessed to produce 
quantum coherent ratchets. 

In this work we consider cold atoms subjected to an 
optical lattice potential that varies smoothly with time, 
instead of being pulsed. We can expect that a driving of 
this kind produces less heating than a kick-type potential, 
and accordingly will preserve the atomic coherence bet- 
ter. We consider a driving potential in which spatial and 
temporal symmetries can be separately controlled, giving 
extreme flexibility for probing and manipulating the sys- 



tem's properties. We firstly show that starting from an 
unbiased initial state, symmetric in both space and time, 
we are able to induce a directed current. Through a per- 
turbative study we find that this ratchet current arises 
from quantum interference between processes which are 
first- and second-order in the driving strength. We then 
provide a simple three-level model to describe its proper- 
ties and find excellent agreement with the exact numer- 
ical simulations. We go on to show that the driving fre- 
quencies at which this current occurs obeys various res- 
onance conditions and discuss the main features of these 



resonances. 



II. ASYMMETRIC DRIVING 

We consider a gas of cold bosonic particles held in a 
toroidal trap. If the lateral dimensions of the torus are 
much smaller than its radius R the system becomes ef- 
fectively one-dimensional. Its low temperature dynamics 
are then well-described by the one-dimensional Gross- 
Pitaevskii equation 



H(t) 



2dx 2 



g\^(x,t)\ 2 + V(x,t) 



(1) 



where the short-range interaction between the atoms is 
described by a mean-field term with strength g. The 
atoms are driven by a time-periodic external potential 
V(x, t) with zero mean, produced by modulating the in- 
tensity of the optical lattice. We note that we measure 
all energies in units of the rotational constant h 2 /mR 2 , 
and set h=l. 

The high level of control available in cold atom exper- 
iments allows us to take the unusual choice of factoring 
the driving potential into separate space and time com- 
ponents 



V(x,t)=KV(x)f(t) 



(2) 



where V(x) gives the spatial dependence of the optical 
lattice potential, and f(t) describes the time-dependence 
of its intensity. The archet ypa l form of a symmetry- 
breaking ratchet potential |12| is V{x) = sin(x) + 
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a sin(2ir, + 4>) , where spatial inversion symmetry is unbro- 
ken for (f> — 7r/2, and is maximally broken for <f> = 0, 7r. 
Accordingly we choose to take 

V(x) — sin(x) + a sin(2ir) , 
f(t) = sin(wi) + /3 sin(2o;i) , (3) 

where the parameters a and /? separately control the spa- 
tial and temporal symmetries of the driving. An exper- 
imental realization of such a driving potential in a cold 
atom system was recently described in Ref. O The 
strength of the driving is denoted by K, and we shall, 
in this paper, restrict ourselves to using small values of 
K for which the system shows a regular response. As K 
is increased the dynamics shows a rich quasiperiodic be- 
havior, and we refer the reader to Ref. [10| for a discussion 
of this and its consequences. 

In Ref. [10| we studied the dynamics under the driv- 
ing ©-(EH when the system was initially prepared in the 
spatially-uniform, time-symmetric state 

V^,0)-(2vr)- 1 /2 . (4) 

This is convenient for experiment as it is the ground state 
of the undriven Hamiltonian, and so can be prepared us- 
ing standard cooling techniques. Clearly the symmetry 
of this state prevents it from carrying a current. We nu- 
merically integrate the wavefunction in time using a split- 
operator method, in each case checking that the time- 
discretization, At, is sufficiently small to produce con- 
verged results. The size of At depended strongly on the 
amplitude of the driving, K, with the surprising result 
that smaller values of K demand a much finer discretiza- 
tion to produce converged results. As an additional ver- 
ification, results were also checked using a fourth-order 
Runge-Kutta technique. To probe the behavior of the 
system we evaluate 

i( t ) ee m)\p*\m) (5) 

as a measure of the current flowing in the ring. We note 
that the expectation value of the momentum operator 
determines the velocity of free propagation in a time- 
of-flight experiment after the optical confinement is re- 
leased. 

Reference [1(3 focused on the case of driving frequency 
u) = 1. It was found that for small K and small or mod- 
erate g the current exhibits large sinusoidal oscillations, 
with a period of ~ 120 driving periods, which clearly av- 
eraged to a non-zero value. As the non-linearity, g, is in- 
creased from zero, the oscillations in current are initially 
enhanced, together with a deformation of the waveform, 
and then become abruptly suppressed above a critical 
interaction strength. This smooth time-periodic behav- 
ior of the current for zero or small interaction strength, 
clearly implies that the driving induces an oscillation be- 
tween the initial state and a single excited state. Exam- 
ining the time evolution in detail, it was found that this 
oscillation occurs chiefly between |0) (the initial state) 



and 1 2), where \l) denotes an eigenstate of the undriven 
Hamiltonian with quantized angular momentum lh. This 
suggests that, for weak driving and weak interactions, we 
can model the dynamics very efficiently by using an ef- 
fective Hamiltonian operating in a reduced Hilbert space 
of a small number of states. The main goal of this paper 
is to study the weak-driving limit in greater detail. The 
resulting perturbative analysis gives an excellent account 
of the system's ratchet dynamics. 

III. FLOQUET STATES: PERTURBATIVE 
STUDY. 

We focus on the behavior of a system governed by the 
the time-dependent Schodinger equation 

i^=Hm (6) 

and we will concentrate on the non-interacting case, 
g = 0, with an initial wave function (j4]). The time- 
periodicity of the Hamiltonian H(t) = H(t + T) implies 
that solutions of the Schrodinger equation are given by 
time-periodic functions known as Floquet states, analo- 
gous to the Bloch wave solutions familiar from studies 
of spatially-periodic systems. Floquet states are of the 
general form 

ipj{x,t) = exp(—iejt)4>j(x,t) , (7) 

where cf>j(x,t) = cf>j(x,t + T) is a time-periodic function 
and the quasienergy Ej is defined modulo 2n/T = uj. 
We also note that our Hamiltonian is of the form 

H(t) = H + V(t), (8) 

where V(t) is given by the driving (UJ)-©. For small K 
we expect that a suitable perturbative study will yield the 
Floquet states from the eigenstates of the unperturbed 
Hamiltonian, which satisfy 

H \l)=e l \l), (9) 

where \l) is an eigenstate of the momentum along the 
ring with space-periodic boundary conditions, with the 
wavefunction 

(x\l) = (2TT)- 1 / 2 exp(ilx) , (10) 

and unperturbed energy e/ = I 2 /2. 

We note first that the time-periodic functions satisfy 
the equation 

[-id t + H{t)]4> j {t)=e j 4> j {t) . (11) 

Since we are interested in the case of space-periodic 
boundary conditions, <f>j(x + 2ir, t) = <fij(x, t) and V(x + 
2ir,t) — V(x,t), we may exploit the combined period- 
icity in x and t (i.e. on the square [0, 2tt] x [0,T]) and 
map the dynamics into a time-independent problem in 
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an effective two-dimensional geometry. Specifically, we 
make the replacement t — > y and rewrite Eq. ([5]) as an 
eigenvalue equation 

[p 2 x /2m+p y + V(x,y)] fo&y) = erfjfay) (12) 

where p y — —id/dy. The generalized Hamiltonian de- 
fined in ([T2"j) is Hermitian. Its only anomaly is that it 
is not bounded from below, but this does not prevent 
us from applying standard tools that do not rely on the 
existence of a minimum energy. So finding the Floquet 
states (J7J) satisfying © amounts to solving the Hermitian 
eigenvalue problem (TT21) . 
We recall that the driving, 

V(x, y) = K[sm(x) + a sin(2x)] [sin(wy) + (3 cos(2wt/)] , 

is small, so a useful starting point is given by the un- 
perturbed Floquet states in the xy representation, which 
satisfy 

(pl/2m + p y )$ m (x,y) = s° lm <lP lm {x,y) , (14) 

where 

cfP lm {f) = (27rT)- 1 / 2 exp( i fc ;m .r) 

= {2TTT)- 1 / 2 exp(ilx - iuj m y) , (15) 

with r = (x,y), and k im = (l,—w m ), where w m = 
2-Km/T — mui, and l,m are integers. The zeroth-order 
approximation to the Floquet quasienergies is given by 

*L = \ ~ "m ■ (16) 

The normalization has been chosen to satisfy orthonor- 
mality, 

f dy P dx [$ m (x 1 y)]'4 m ,(x,y)=8 lv 6 mm , . (17) 
Jo Jo 

In the following we focus on those values of I and m 
which satisfy the resonance condition 

p 

u) m = — , i.e. eL = , (18) 

where we expect to find the highest values of the ratchet 
current. 

The periodic driving (|13j) has a finite number of Fourier 
components, so that in the expansion 

V(f)=Y,Vjtexp(i3-t f ) (19) 
S 

only a handful of reciprocal lattice vectors g satisfy 
Vg 7^ 0. Specifically, the driving Fourier component is 
nonzero for the 16 values of g represented by small red 
dots in any of the four figures shown in Fig. Q] In each 
of these figures, the blue circles indicate the ki m vectors 



that satisfy the resonance condition (fT8|h which can also 
be written as 

to = l 2 /2u , (20) 

and whose continuous version is represented by the green 
parabola. 

In this peculiar xy space all resonant states are "degen- 
erate" with energy e° m = 0. These states are connected 
by the driving (p~3|) . (fT9|) . By construction, our initial 
state corresponds to the unperturbed Floquet state |00), 
where in general \lm) is the state characterized by ki m , 
its wave function being 

(x,y\lm)=$ m (x,y) . (21) 

This state |00) is represented by the filled blue circle in 
each of the four graphs of Fig. [1] Thus we expect the 
driving to mix |00) with the other resonant Floquet states 
satisfying = 0. 

Let us focus first on the case u = 1 (upper-left graph 
of Fig. [I}. Starting from the state 

^ (x s y) = (27 I T)-V2 1 (22) 

we expect the system to evolve from 1 00) towards 1 22) and 
1 22), where I = —I. In general, under the effect of driv- 
ing the state |00) will also mix with higher- lying unper- 
turbed Floquet states, but we may expect that mixing to 
be weaker as it involves higher powers in the Fourier com- 
ponents Vg, all of which are small because K is assumed 
to be small. Thus a truncated Hilbert space spanned 
only by the Floquet states {|00), |22), |22)} may suffice 
to describe the dynamics under weak driving. The next 
section is devoted to the formulation of the three-level 
model. 



IV. THREE-LEVEL MODEL: EFFECTIVE 
MATRIX ELEMENTS. 

Inspection of the upper-left Fig. Q] strongly suggests 
that for uj = 1, the mixing between |00) and e.g. the 
state 1 22) must be dominated by the interference between 
the first-order process involving g = (2, 2) once, and the 
the second-order process involving g = (1,1) twice (with 
g given in the units of the array there shown), while vis- 
iting the state |11) virtually. In this section we calculate 
the effective matrix elements which account for this first- 
and second-order mixing in an effective time-independent 
problem defined in the xy space. 

First we note that, in this xy representation, we deal 
with a formally time-independent problem defined by the 
Hamiltonian 

H = H Q + V , (23) 

where Hq is given in (IT41 and V by (fl"3"|) . The full Green 
function G(z) = (z — H) -1 can be related to its un- 
perturbed counterpart G°(z) = (z — Hq)^ 1 through the 
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FIG. 1: Schematic representation of the driving-induced transitions for driving frequencies uj = 1/4, 1/2, 1,2; see main text 
for the meaning of the various symbols. In each case the direct, first-order (solid arrows) and indirect, second-order (dashed 
arrows) transition amplitudes interfere to produce a ratchet current. For clarity, in the u) = 1/2 figure the first dashed arrow 
of the various second-order processes has been omitted. 



Dyson equation \l£ 



G = G Q + G°VG 



(24) 



One may also describe the effect of the perturbation V 
in terms of the T-matrix (not to be confused with the 
time-period) satisfying 



or, equivalently, 



T = V + VG°T 



G = G° + G°TG°. 



(25) 



(26) 



Thus the effective dynamics up to second order in V may 
be described in terms of an effective T-matrix approxi- 
mated as 



T{z) ~ V + VG°{z)V 



(27) 



Finding the matrix elements of this approximate, second- 
order T-matrix is equivalent to finding those of an effec- 
tive V which, treated to first order, yields the correct 



second-order dynamics. So the next goal is to compute 
the matrix elements (lm\T\l'm'). Applying the closure 
relation to (l27l) we obtain 



(j\T(z)\j') = (j\V\j')+J2 



(j\V\j")(j"\V\f) 



(28) 



where j is a short-hand notation for the quantum num- 
bers Im. 

For the case lj = 1 we restrict our analysis to 
the matrix elements between the three resonant states 
{|00), 1 22), 1 22)}. These three states have an unperturbed 
energy e° = 0, so in (|2"5]l we focus on the shell z = 0. In- 
spection of Fig. Q] clearly shows that the intermediate 
state connecting |00) and |22) is |11), while |11) is 
the intermediary between |00) and |22). Thus, for exam- 
ple, 



(00|T|22) = (00|V|22)- 



(00|V|11)(11|V|22) 



= v 22 - 



r 2 



5 



K 

T 



-a/3- 



~2 



(29) 



where V; m stands for the Fourier component Vg, with 
g = (/, m) given in the units of Fig. [T] 
Similarly, 



with 



(00|T|22) 



K o K 



(30) 



(31) 



Restricting to these two well-behaved second-order ma- 
trix elements, we construct a three-level model defined by 
the Hamiltonian matrix 



Ho 



o r 4 



o 



r + o r_ 
o r_ o 



(32) 



spanning the space {|22), |00), |22)} or, for brevity, 
{|2), |0), 1 2)} (using this ordering). In the subspace 
{|2), 1 2)} it is always possible to introduce a rotation such 
that one state is decoupled from |0). Specifically, if we 
define the orthonormal states 



1 



\b) 



1 



(r+|2) + r_|2» (33) 
(r_|2)-r+|2)) , (34) 



we find that (0|T|6) = 0, so that p2|) transforms into 



o r o 
roo 
ooo 



where 
T = (0\T\a) = 



r 2 _ = 



K 



(35) 



V^ 2 +4a 2 /3 2 . (36) 



Therefore, for each particular driving the general three- 
level model can be truncated to an effective two-level 
problem, as found numerically in Ref. [ToL 



V. AVERAGES OF THE CURRENT 

Returning to the real time (xi) picture, it is clear that 
if the system is initially prepared in a given state ip(x, 0), 
its subsequent evolution for t > 0, when the driving is on 
[we assume ([2]) is multiplied by a step function 8(t)], can 
be economically written as an expansion in the complete 
basis of the Floquet states, of general form ( [7| 

i>(x,t) = ^2cjexp(-iejt)(l)j(x,t) , (37) 



= / ^) j (x ) 0)*tp(x,0)dx , (38) 



In this sense the Floquet states represent a generalization 
of the standard energy eigenstates of static Hamiltonians 
to the case of time-periodic systems. 

Within our perturbative approach, we have seen that 
the three-level picture reduces in each case to a two-state 
(|0) and \a)) problem, of Hamiltonian (|35|) . Thus, if the 
initial state is | V(0)> = |0), or tp(x,0) = (2-kT)- 1 / 2 , the 
system will undergo a simple oscillation of the Rabi type 
between the unperturbed Floquet states |0) and \a). The 
frequency of those oscillations will be 2r. Oscillations 
of exactly this type were observed in the numerical in- 
vestigation of this model in Ref. [l^ and in the recent 
experimental investigation (l3| . 

As well as evaluating the time-dependent current, it 
is also convenient to calculate the time- averaged current. 
This is particularly important to verify the existence of a 
long-lasting ratchet effect, which requires that the time- 
averaged current remains non-zero as the observation pe- 
riod tends to infinity. Such time-averages can be formed 
in two distinct ways. The first is the stroboscopic average, 
in which the current is evaluated only at discrete times 
t n = to + nT, where t G (0, T}. It can be defined as 



Is(t ,N) 



1 N 



N 



nT) 



(39) 



n=0 



This is frequently the scheme of measurement most con- 
venient for experiment. One example is when the driving 
potential V(x, t) is obtained by periodically accelerating 
and decelerating the optical lattice, thereby producing 
an inertial force in the rest frame of the lattice. Mea- 
surements of the current, however, are made in the rest 
frame of the laboratory, and so it is convenient to make 
measurements at times when the laboratory and lattice 
rest frames coincide, which occurs stroboscopically. The 
other averaging scheme is a continuous time average, 



Ur) = T 



I(t)dt 



(40) 



When the period of the driving is much shorter than the 
response of the system, these averages will coincide at 
long times. However, for lower driving frequencies it is 
possible that the two time-averages will yield different 
results, as the stroboscopic sampling will only capture a 
subset of the system's dynamics. 

At long times (JV, r -> oo), we can use ©, f3"7]>. (|3"9"1) . 
and (|40p to write these time-averaged currents as 

T s (t ) = 2l c il 2 ^i(*o)IP*l^(*o)) ' ( 41 ) 



Ic = ~ / dt I s (t ) 
1 Jo 



(42) 



In (j4T]) we have assumed that e 3 7^ e.y for all j 7^ j'. 

Once the driving has been switched on and Rabi oscil- 
lations have started, the system will spend, on average, 
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half its time in |0) and \a) . Thus the analytical prediction 
for the ratchet current is simply (7 = J c ) 



I = h/2 



(43) 



where I a = (a\p x \a) is the current carried by \a). 

We first consider the case of the main resonance to = 1 . 
In this case the Rabi oscillation occurs between |0) and 
| a) given by ([33]) , with frequency 



9., 



2T 



K 



y/K 2 + Aa 2 /3 2 



(44) 



We note that, in our reduced units, the states |2) and 
|2), carry currents 2 and -2, respectively. As a result, the 
ratchet current can be shown to be 



I = 



4Ka/3 



K 2 +4a 2 f3 2 



(45) 



From (|45l) the beautiful picture emerges of a coherent 
ratchet current stemming from the interference between 
first- and second-order processes creating an imbalance 
between the matrix elements coupling the initial, time- 
symmetric state |0) to the current-carrying states |2) and 
|2). Most importantly, the ratchet current exists only if 
both a and (3 are nonzero, i.e. if the driving is both space- 
and time-asymmetric, in agreement with the symmetry 
analysis provided in [15(. We can further observe that 
the ratchet current is a unique function of the product 
a/3, and for this driving frequency the ratchet current is 
maximized for a/3 = K/2, in excellent agreement with 
the experimental observation [l3j . 

Figure [2] shows the extremely good agreement between 
the analytical prediction of the three-level model and the 
full numerical simulation for the time-dependent ratchet 
current in the weak coupling case (K = 0.1). In partic- 
ular, this means that the analytical predictions (|4"4"]) and 
(1451) for the Rabi frequency and the long-time ratchet cur- 
rent become essentially exact in the weak-driving limit. 
As expected, the agreement with the perturbative ana- 
lytical calculation worsens for higher K. 




50 
time [T] 



FIG. 2: Comparison of the time-dependent current predicted 
by the effective three-level model and the exact numerical 
results, for ui — 1 and asymmetry parameters a = P = 0.2. 
The three-level model predicts Rabi oscillations which fit the 
exact results extremely well for weak driving (K = 0.1). For 
strong driving (K = 0.5) the exact results show additional 
small, high-frequency oscillations, but the main behavior is 
still reasonably well-described by the effective model. 



an enlargement of the small-scale structure in the current 
response, which indicates the existence of further fami- 
lies of sub-resonant peaks, of much smaller magnitude 
than the four primary peaks, occurring at commensu- 
rate fractions of the driving frequencies. This perhaps 
indicates the role of higher order interference processes, 
which could in principle be described by a suitable gen- 
eralization of our procedure. 



We recall that our analysis predicts, in addition to 
the main resonance at ui = 1, additional resonances for 
a; = 1/4, 1/2, and 2. To verify this prediction we show 
in Fig. [3] the time-averaged current obtained for a fixed 
driving strength of K = 0.05 as the frequency is var- 
ied over a wide range. We can first note that we indeed 
see peaks at the four resonant frequencies indicated by 
our model. A similar set of resonance peaks was ob- 
served in the experimental investigation of this ratchet 
(l3j . The ui = 1 peak is considerably larger than the 
others, and it is interesting to note that the response for 
uj = 1/2 is of opposite sign to the other peaks. This peak 
is also unusual in its sensitivity to the averaging proce- 
dure used; the stroboscopic result is larger and broader 
than the continuous time-average. In the inset we show 



In the same way as for the u = 1 resonance, we can ob- 
tain analytical results for the ratchet currents produced 
by the other main resonances, governed by first and sec- 
ond order transitions as indicated in Fig. [1] These results 
are given in Table HI In Fig. 0] we compare the analyti- 
cal and numerical predictions for the continuously time- 
averaged ratchet current for the four main resonances. 
The agreement is excellent. We note in particular the lin- 
ear behavior I oc K for small K, which is also predicted 
analytically, and the existence of a maximum current for 
w = l, occurring for a driving strength of K = 0.08. 
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FIG. 3: Time-averaged current, averaged over 4000 driving 
periods for a weakly-driven system (K — 0.05, asymmetry 
parameters a = ft = 0.2), plotted as a function of the driv- 
ing frequency. Four peaks appear, at driving frequencies of 
oj — 0.25,0.5, 1,2 in agreement with the resonant condition 
(|18|) . Inset: Magnified view, showing the existence of sub- 
resonances at commensurate fractions of the resonant frequen- 
cies. We also note the reduction of the to — 0.5 resonance 
when the time-average is evaluated continuously instead of 
stroboscopically. 




driving strength K 



FIG. 4: Comparison of the analytical results with 
numerically-exact data for the four principal resonances, ui = 
0.25, 0.5, 1,2 8 (see Table I). The values for uj — i are actu- 
ally negative, and so for convenience we plot their absolute 
value. For all the curves the asymmetry parameters were 
set to be a — /3 — 0.2. We can see that for weak driving 
strengths the agreement is excellent; the effective three-level 
model produces quantitatively accurate results. For higher 
driving strengths the model diverges from the exact results, 
as expected for a perturbative result. 

As well as the amplitude of the ratchet current, the ef- 
fective model also provides predictions for the period of 
the Rabi oscillations. In Fig. [5]we compare these predic- 
tions [see Eq. with the numerically exact results, 
and again see excellent, quantitative agreement over a 
wide range of driving amplitudes. 



Finally, in Fig. [5] we show the stroboscopically- 
averaged current as a function of the initial sampling 
point time to. For to = 0.25, 1, and 2 the stroboscopic 
quantity shows only a very weak dependence on ty. The 
uj = 0.5 resonance, however, does display an important 



0.001- 

om ' o'.i 

K 

FIG. 5: Rabi frequencies, Q,r, predicted by the effective three- 
level model [see Eq. (|44p and Table U are shown with solid 
lines, to compare with data extracted from numerically exact 
simulations. As in Fig. fj] we see excellent agreement between 
analytics and numerics. 




time [T] 

FIG. 6: The stroboscopic average is evaluated at discrete 
times t n — to + nT. Here we show the dependence of this 
time average on the initial time sampling point to. For 
uj — 0.25, 1, 2 the variation is weak, but for uj = 0.5 the 
dependence is strong, and contributes to the low value of 
the continuous current average for this value of the driving 
frequency. The number of driving periods used to compute 
the stroboscopic average is, from top to bottom, 400, 400, 
600, 300. This amounts to about ten complete Rabi periods 
(2-k/Q.r) in the first two cases, and one Rabi period in the 
last two. 



dependence on to, and indeed changes sign as to varies 
over the range < to < T. This causes the continuously- 
averaged current to be significantly smaller than a strobo- 
scopic estimate for this driving frequency [see Eq. (|42|) ]. 
as can be seen in Fig. [3] 
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LU 0_R I 

1 K h , / 1 a\2~\ 1 / 2 lOKa/3 

1 172° Y K + ( 2a ^ \ X* + (2aV 
o K \l\ tsq\2 . „21 1 / 2 70Ka^ 

2 75 LjjfjjP) + Q J (35a)*+(W 

TABLE I: Analytic results for the Rabi frequencies and the 
continuously time-averaged ratchet currents for the four res- 
onant frequencies considered. 



VI. CONCLUSIONS 

We have studied an unusual form of flashing ratchet, in 
which the spatial and time symmetries can be controlled 
independently. Using a novel form of perturbation the- 
ory, we find that we are able to describe the weak-driving 
regime of this system to a surprisingly high degree of 



accuracy by using a three-level effective model, which 
later reduces, in each particular case, to a simple two- 
level model. This provides an analytical underpinning to 
the phenomenological two-level model introduced in Ref. 
[Tol to describe this ratchet system. The ratchet current 
arises from the interference between first- and second- 
order driving-induced processes, and so is a purely quan- 
tum coherent effect, not involving dissipation. It should 
be noted that we have neglected the effect of the non- 
linear interaction g. In Ref. 10 it was shown that this 
has the effect of damping the Rabi oscillations, eventu- 
ally producing a "self-trapped" state. Introducing the 
non-linearity in a consistent way to this model remains a 
interesting topic for future work. 
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